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We discuss the angular resolution obtained for events registered with the surface detector alone and for hybrid 
events, i.e., those observed simultaneously by both the surface and fluorescence detectors. The angular accuracy 
of the surface detector is directly extracted from the data itself and on an event by event basis, and is given as 
a function of the number of stations triggered by the event and of the zenith angle of the shower. We compare 
the angular resolution of the surface detector obtained from hybrid events with the one obtained from the surface 
detector alone. 



1. Introduction 

To fulfill the goal of the Pierre Auger Obser- 
vatory of determining the origin of the ultra high 
energy cosmic rays, their arrival direction must be 
obtained with optimal accuracy. This includes a 
precise measure of the angular resolution (AR), 
which we define as being the angular radius that 
contains 68% of the showers coming from a given 
point source. 

The Pierre Auger Observatory consists of 
two independent components: the surface de- 
tector (SD), which is comprised of 1600 water 
Cherenkov detectors distributed in a triangular 
grid with a separation of 1.5 km and cover- 
ing an area of 3000 km 2 ; and the fluorescence 
detector (FD), which is formed by 24 fluores- 
cence telescopes located in 4 buildings overview- 
ing the surface detector pQ. Hybrid events, i.e., 
events observed simultaneously by both compo- 
nents, have smaller reconstruction uncertainties 
than the events observed with only the surface 
component. On the other hand, the latter have 
much higher statistics than the former. 

The angular resolution for hybrid events is de- 
termined from simulations, by computing the an- 
gle between the simulated event and its recon- 
structed direction (section [2]). The angular reso- 
lution for the SD is determined, on an event by 
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event basis, from the zenith (9) and azimuth (0) 
uncertainties obtained from the geometrical re- 
construction (section [3|) , using the relation [2|3j : 

F(rj) = 1/2 (V[0\ +sin 2 (6») V [(/>]) 

where r\ is the angle between the reconstructed 
shower direction and the true one, and V[9] and 
V[4>] are the variances of 9 and 4> respectively. If 
9 and <fi/ sin(6) have Gaussian distributions with 
variance cr 2 , then F(rf) — a 2 and r] has a distribu- 
tion proportional to e~ r ' l 2a d(cos(rj)) d(f>. Then, 
according to our definition, the angular resolution 
is AR = 1.5 y/F(rj). 

The arrival direction of a SD event is deter- 
mined by fitting the arrival time of the first par- 
ticle in each station according to a shower front 
model. The accuracy achieved in the arrival di- 
rection depends on the clock precision of the 
detector and on the fluctuations in the arrival 
time of the first particle. The timing uncertainty 
is directly modeled from the data in each sta- 
tion [4j (section l3~Tj) . The model is adjusted us- 
ing pairs of adjacent stations located in the sur- 
face array. Then, it is validated by studying the 
X 2 probability distribution for the geometrical re- 
construction and by comparing two independent 
reconstructions of the same event (section 13 .2[) . 
The angular resolution can therefore be estimated 
for the SD reconstruction on an event by event ba- 
sis ( section [3. 3p . Using the hybrid events, we are 
able to extract, for a subset of events, the angular 
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resolution of the surface detector and compare it 
with the one obtained on an event by event basis 
(section 0} . 

2. Angular resolution of hybrid events 

The angular resolution for hybrid events is de- 
termined from simulations, by computing the an- 
gle (77) between the injected shower axis and the 
reconstructed one. We simulated 5678 proton 
showers with the Corsika 6.616 code [5] using 
QGSJetII-03 [6J and FLUKA [7\ codes for high 
and low energy hadronic interaction models re- 
spectively. The showers were simulated up to a 
zenith angle (9) of 65°, with a distribution of 
cos(#) sin(0), an uniform distribution in the az- 
imuth angle, and for various energies from 0.1 
EeV to 10 EeV in steps of 0.25 in logarithmic 
scale. For each shower we generated uniformly a 
core position in a slice of 60° in front of one of the 
fluorescence telescopes with a maximum distance 
from the telescope increasing with the energy ac- 
cording to the trigger efficiency [8] . Once the hy- 
brid simulation was performed, we reconstructed 
the simulated events with the same reconstruc- 
tion chain used for data. 

In figure [T] we show the angular resolution, de- 
termined as the value where the cumulative distri- 
bution function of rj{9) reaches 0.68, as a function 
of the energy. As can be observed, the angular 
resolution improves for larger energy showers. It 
is worthwhile to remark that this behavior origi- 
nates from the convolution of the dependence of 
the angular resolution with the different geomet- 
rical parameters. Also, it is important to check 
these results using real data. For example, using 
stereo events, which have fluorescence telescopes 
with triggered pixels in more than one building. 

3. Angular resolution of the surface detec- 
tor 

3.1. Time variance model 

The angular accuracy of the surface detector 
events is driven by the precision in the measure- 
ment of the arrival time of the shower front (T s ) 
at each station [3]. The particle arrival time in 
the shower front can be described as a Poisson 
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Figure 1. Angular resolution of hybrid events de- 
termined with simulations as a function of the 
energy. 

process over a time interval T. The arrival time 
of the first particle (Ti) is used as the estimator 
for the shower front arrival and its distribution 
function is given by: 
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with t = T/n, where n is the number of particle^) 
measured during the time T. Then, the variance 
of Ti is t 2 , but since we estimate the parameter 
T from the data itself, it is modified to [I]: 
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The variance of T s in the SD stations is given by 
the sum of the detector clock precision (b 2 ) and 
of the variance of T%. It then becomes: 



V[T a \ 
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where T50 is the time interval that contains 
the first 50% of the total signal as measured 
by the photomultiplier FADC (flash analog-to- 
digital converters) traces. The parameter a is a 



2 To calculate the number of particles n, we assumed that 
the muons are the particles that contribute the most to 
the time measurement and in average cross the detector 
with the direction given by the zenith angle of the shower. 
The number of particles is calculated as the ratio between 
the total signal in the stations and the track length of the 
particles. 
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scale factor, containing all the assumptions con- 
sidered in the model and the treatment done to 
the FADC traces. The parameter b should be 
given by the GPS clock accuracy (about 10 ns) 
and the FADC trace resolution (25/VT2 ns), that 
is b ~ 12 ns. Both a and b are determined from 
the data. 

A special sub-array of pairs of water Cherenkov 
detectors has been deployed as a part of the sur- 
face array. These are adjacent detectors located 
~11 m apart, and therefore arc sampling the same 
region of the shower front. To determine the pa- 
rameters a and b we used all the events that pass 
our selection criteria [9] from April 2004 to the 
end of August 2008 with at least one pair in the 
event. There is total of 46416 events, which are 
used to fit these two parameters yielding: 

a = 0.60 ±0.01, 
b = (14.1 ± 0.2) ns, 

with a x 2 /ndof — 0.9992. More details about the 
time variance model and the fitting procedure can 
be found in reference [4|. 

3.2. Validation of the time variance model 

In this section, we show that the model cor- 
rectly reproduces the uncertainties of the arrival 
time of the first particle in the stations. 

We define the time difference AT = dT x - dT 2 , 
where dTj (dT^) is the time residual of the first 
(second) twin station to the fitted shower front. 
If the time variance model describes correctly 
the measurement uncertainties, the distribution 
of AT/-y/V[AT], where V[AT] = V[T{\ + V[T 2 ], 
should have unit variance. In figure[2]we show the 
RMS of the distribution of AT/ ^V[AT] for the 
adjacent detectors as a function of cos(#) (top), 
the average signal (middle), and the distance of 
the paired detectors to the core position (bot- 
tom). In all the cases, the RMS is almost con- 
stant and close to unity, which shows that the 
time variance model adequately reproduces the 
experimental data. It is worthwhile to notice that 
the time variance model does not explicitly de- 
pend on the distance of the station to the shower 
core. Despite this, the result shown in the bottom 
panel is satisfactory, with only a small tendency 
to overestimate the variance for detectors very 
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Figure 2. T he RMS of the distribution of 
AT/ y/V[AT], as a function of the shower zenith 
angle (top), the average signal in the adjacent de- 
tectors (middle) , and their distance to the shower 
core (bottom). 

near to the shower core. 

We also studied the distribution of the x 2 prob- 
ability of the geometrical reconstruction. In fig- 
ure [3] (top) we show the distribution for the 
308234 events passing our quality cuts [9] with 4 
or more stations. This distribution is almost flat 
as it should be in the ideal case, despite the small 
peak at low values. This small peak could be due, 
for example, to some stations in the events that 
have low signal, causing a large x 2 ■ But as is seen 
in figure^ (bottom), the peak disappears for large 
multiplicity events (for 5 or more stations). For 
both cases, we only plot \ 2 probabilities larger 
than 1% to avoid the large peak at zero corre- 
sponding to badly reconstructed events (~ 5%). 
Also, we study the x 2 probability distribution for 
two zenith angle ranges, as shown in figure [U The 
X 2 probability distribution is flat for both large 
and small zenith angles, which means that the 
model works for all angles without compensating 
one set from the other. This distribution shows 
that the variance model properly reproduces the 
uncertainties of the arrival time of the particles 
in the stations and allows us to determine the 
angular resolution from the uncertainties in the 
reconstruction data. 

We are able to validate independently the time 
variance model by using the redundant informa- 
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Figure 3. The x 2 probability distribution for all 
events (top) and for 5 or more stations (bottom). 
See text for details. 
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Figure 5. Distribution of the difference of zenith 
(left) and azimuth (right) angles from each in- 
dependent reconstruction divided by the square 
root of the sum of its variances for the case of 6 
or more stations in the event. 
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Figure 4. The x 2 probability distribution for 
events with 5 or more stations and with zenith 
angle smaller (top) and higher (bottom) than 55°. 



tion given by the special sub-array of adjacent 
detectors. To perform this analysis, we select 
those events with at least three pairs and with the 
shower core inside the region defined by these sta- 
tions, to guarantee a good reconstruction. Then 
the reconstruction is performed twice, each time 
using the information of one of the tanks in each 
pair. This provides two quasi-independent esti- 
mates of the arrival direction of the same shower. 
In the left (right) panel of figure [5] we show the 
distribution of the difference of the zenith (az- 
imuth) angles from each reconstruction divided 



by the square root of the sum of its variances for 
the case of 6 or more stations which corresponds 
to events with energy higher than 10 EeV. We fit 
the distribution with a Gaussian and we obtain 
a sigma value compatible with 1. Despite of the 
low statistics, this result shows that the uncer- 
tainties in the determination of the arrival time 
of the first particle in the stations are well repre- 
sented. The same results are obtained for other 
multiplicities, except for the 3-fold case, where 
the sigma is about 0.75 which indicates that we 
are slightly overestimating the uncertainties. 

3.3. Angular resolution 

Considering the quality of the time variance 
model for the measurement uncertainties, we can 
calculate directly the angular resolution on an 
event by event basis out of the minimization pro- 
cedure. In figure [6l we show the angular resolu- 
tion as a function of the zenith angle for various 
station multiplicities. 

The angular resolution is about 2.2° in the 
worst case of vertical showers with only 3 sta- 
tions hit. This value improves significantly for 4 
and 5 stations. For 6 or more stations, which cor- 
responds to events with energies above 10 EeV, 
the angular resolution is in all cases better than 
1°. Above 60°, the event multiplicity increases 
rapidly with zenith angle, and only a few low en- 
ergy events trigger only 3 stations, thus the ac- 
curacy decreases. Also, a hump appears around 
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Figure 6. Angular resolution as a function of 
the zenith angle {&) and for events with 3 sta- 
tions (circles), 4 stations (squares), 5 stations 
(up-triangles), and 6 or more stations (down- 
triangles). 

40°, more visible in the 3-fold case. This is due 
to the contribution of the uncertainties on the 
core position for events with lower energy (less 
than few EeV) . In figure we show the angu- 
lar resolution as a function of the zenith angle 
for events with an energy above 3 EeV, which is 
the energy corresponding to a trigger efficiency 
greater than 99% [TO]. In spite of the degrada- 
tion of the statistics for 3-fold events, the hump 
disappears and the angular resolution gets better. 
For high multiplicity events, as it is expected, the 
angular resolution is not affected by the cut in 
energy. We want to remark that all uncertainties 
quoted in figures d] and [7] are statistical only. We 
did not, at this stage, investigate possible biases 
or systematics in the determination of the arrival 
direction angles. 

4. Angular resolution of the surface detec- 
tor using hybrid events 

Hybrid events that trigger 3 or more stations 
can be reconstructed using both the hybrid and 
the surface detector alone modes, giving two in- 
dependent estimates of the geometry. The com- 
parison of these estimates is therefore an addi- 
tional independent check of the accuracy on the 
determination of the arrival direction of the cos- 
mic rays. Therefore, we compute the angle (77) of 
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Figure 7. Angular resolution as a function of the 
zenith angle (9) for events with an energy above 
3 EeV, and for various station multiplicities. 

those two estimates for different mutiplicities and 
zenith angle ranges. Then, the angular resolution 
of the surface detector can be obtained as: 

ar sd = Jar* ARi yb 

where AR V is the value of rj for which the cumula- 
tive distribution function of r](0) reaches 0.68 and 
ARhyb is the angular resolution of hybrid events 
obtained from Monte Carlo simulations (shown in 
figure [JJ . 

In figure [H] we show the angular resolution for 
the surface detector (ARsd ) using the hybrid re- 
construction as a reference, as a function of the 
zenith angle (6) for different numbers of stations 
in the event. This plot can be directly compared 
with the one in figure [5J despite the differences in 
the statistics. The values obtained in figure [S] are 
slightly higher than the ones obtained from the 
direct determination of the angular resolution on 
an event by event basis. This could be due to sys- 
tematic uncertainties in the reconstructions. For 
the case of 3-fold events, this slight difference is 
not present and this could be due to the fact that 
the uncertainties seem to be overestimated (sec- 
tion [32]). 

5. Summary 

We have determined the angular resolution for 
the hybrid events using Monte Carlo simulations 
and we found that for energies above 3 EeV it 
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Figure 8. Angular resolution for the surface de- 
tector events using the hybrid reconstruction as 
a reference, as a function of the zenith angle (8) 
and for different numbers of stations in the event. 
See text for more details. 

is about 0.5°. For the surface detector, we de- 
termined the angular resolution from the data. 
This could be done thanks to the time variance 
model developed to describe the measurement un- 
certainties of the time arrival of the first parti- 
cle in the surface stations. Using this model we 
obtaines an optimal determination of the shower 
arrival direction and are able to extract the an- 
gular resolution of the surface detector on an 
event by event basis. We found the angular res- 
olution to be better than 2.2° for 3-fold events 
(E < 4 EeV), about 1.5° for 4-fold and 5-fold 
events (3 EeV< E < 10 EeV) and better than 1° 
for higher multiplicity events (E > 10 EeV). 
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